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ABSTRACT 

Current measurements indicate that the deuterium abundance in diffuse interstellar gas varies spatially 
by a factor of ^4 among sightlines extending beyond the Local Bubble. One plausible explanation for 
the scatter is the variable depletion of D onto dust grains. To test this scenario, we have obtained high 
signal-to- noise, high resolution profiles of the refractory ion Ti II along seven Galactic sightlines with D/H 
ranging from 0.65 to 2.1 x 1CP 5 . These measurements, acquired with the recently upgraded Keck/HIRES 
spectrometer, indicate a correlation between Ti/H and D/H at the > 95% c.l. Therefore, our observations 
support the interpretation that D/H scatter is associated with differential depletion. We note, however, 
that Ti/H values taken from the literature do not uniformly show the correlation. Finally, we identify 
significant component-to-componcnt variations in the depletion levels among individual sightlines and 
discuss complications arising from this behavior. 
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In standard big bang nucleosynthesis (BBN), deuterium 1 At r nt it ms th t th D/H ratio has 

is the most sensitive baryometer of the light elements (Schramm cer am ieS ''' , , P re f? n TTT a 6 ' ra 10 as 

intrinsic scatter within the Milky Way. 



& Turner 1998). Measurements of its primordial value 
in quasar absorption line systems have placed tight con- 
straints on the baryonic mass density fib (Buries & Tytler 
1998; O'Meara et al. 2001; Kirkman et al. 2003). The 
inferred fib value is in impressive agreement with the value 
derived from CMB measurements of the WMAP (e.g. Spergel 
et al. 2003) and other microwave experiments lending con- 
fidence to both the deuterium measurements and BBN 
theory In addition to its cosmological significance, deu- 
terium is an important tracer of chemical evolution. Deu- 
terium is astrated within stellar cores, and there are no 
known means of producing significant amounts of D other 
than BBN. The evolution of D/H, therefore, tracks the 
global history of the gas astration within a galaxy (e.g 
Chiappini, Renda, & Matteucci 2002). 

It is in this context that measurements of D/H within 
the Milky Way have impact: (1) the upper bound to the 
Galactic D/H value sets a lower limit to the primordial 
D/H value; and (2) a comparison of Galactic D/H with the 
primordial D/H value describes the chemical evolution his- 
tory of the Galaxy. Since the first measurements of Galac- 
tic D/H with Copernicus (Rogerson & York 1973), sur- 
veys for Galactic D /H have been pursued on each succeed- 
ing UV observatory bearing a high resolution spectrometer 
(e.g. Linsky et al. 1995; Jenkins et al. 1999; Hoopes et al. 
2003; Wood et al. 2004). The most remarkable conclu- 
sion of these efforts is that beyond the Local Bubble the 
D/H values have significant scatter; D/H ranges from 5 to 
20 x 10~ 6 . This dispersion is larger than the estimated 
statistical error and is unlikely to be systematic error as- 
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Recently, observations with the Interstellar Medium Ab- 
sorption Profile Spectrograph (IMAPS) and the FUSE Ob- 
servatory have extended D/H measurements to large dis- 
tances from the Sun and correspondingly higher iV(HI) 
values. Although only a handful of measurements ex- 
ist to date, Wood et al. (2004) and others have noted 
that the few sightlines that probe the greatest distances 
(> 500pc) have a central value and dispersion that are sig- 
nificantly smaller than sightlines at intermediate distances 
(20 pc < d < 500 pc). One possible interpretation of these 
trends is that D is significantly depleted onto dust grains 
(e.g. Jura 1982). Draine (2004a, b) has examined the 
principal mechanisms of D adsorption and grain destruc- 
tion and argues that it is at least plausible that D would be 
preferentially depleted. He proposed testing the depletion 
hypothesis by comparing the D/H values with abundances 
of refractory elements like Fe, Ni, and Ti (i.e. species that 
are highly prone to depletion onto grains) along the same 
sightlines. 

We have initiated a program to obtain high signal-to- 
noise (S/N), high resolution observations of Ti II profiles 
for sightlines with accurate Galactic D/H measurements. 
Previous surveys of Ti have demonstrated that it is highly 
refractory, presumably because of its large condensation 
temperature (e.g. Stokes 1978; Lipman & Pettini 1995). 
Therefore, a measurement of Ti/H in the ISM assesses the 
depletion level along that sightlinc. Of additional impor- 
tance, Ti + is the dominant ion in H I regions. Unlike 
Na° and Ca + the ionization potential of Ti + is rs 1 Ryd 
and Ti + is predominantly shielded from ionizing photons. 
Furthermore, Ti + has a charge exchange reaction rate with 
hydrogen that is significantly greater than many other ions 
(e.g. Fe + ,Si + ; Kingdon & Ferland 1996). We expect is less 
sensitive to photoionization effects and Ti II should trace 
the velocity profiles of H I and D I gas. Therefore, high 
resolution Ti II profiles are likely to be better suited than, 
e.g., Fe II, for constraining the fits of D I and assessing 
the likelihood of deuterium line saturation in high AT (HI) 
sightlines. 
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Table 1 

SUMMARY OF OBSERVATIONS 



Target 


eOri 


5 Oil 


HD 191877 


HD 195965 


BD+28°4211 


tOri 


FcigcllO 


B 


1.71 


2.01 


6.22 


6.60 


13.01 


2.53 


13.01 


Obs (UT) 


08Sep2004 


08Sep2004 


08Sep2004 


08Sep2004 


06Oct2004 


08Sep2004 


09Sep2004 


tezp (s) 


9 


13 


1800 


1200 


1800 


9 


1800 


c /Ma 

o/JN" 


OTA 

370 


355 


OCA 

350 


250 


145 


AO A 

280 


120 


log[iV(HI)/cni- 2 ] 
D/H h 


20.40 ±0.08 


20.19 ±0.03 


21.05 ±0.05 


20.95 ±0.03 


19.85 ±0.02 


20.16 ±0.10 


20.14 ±0.07 


0.646 ±0.38 


0.736 ±0.12 


0.776 ±0.20 


0.851 ±0.17 


1.380 ±0.10 


1.413 ±0.88 


2.138 ±0.43 


Ref 


1 


2 


3 


3 


4 


1 


5 


log[7V(TiII 3230)] 


< 11.55 


< 11.46 


12.26 ±0.03 


11.97 ±0.05 


< 11.57 


< 11.67 


< 11.77 


log[7V(TiII 3242)] 


11.50 ±0.05 


11.28 ±0.06 


12.25 ±0.02 


12.06 ±0.02 


< 10.99 


11.20 ±0.07 


11.69 ±0.05 


log[iV(TiII 3384)] 


11.37 ±0.04 


11.09 ±0.06 


12.23 ±0.02 


12.00 ±0.02 


11.08 ±0.08 


11.34 ±0.03 


11.54 ±0.04 


log[iV(TiII)] c 


11.40 ±0.03 


11.15 ±0.04 


12.24 ±0.02 


12.02 ±0.02 


11.08 ±0.08 


11.31 ±0.03 


11.59 ±0.03 


Iog(Ti/H) 


-9.00 ±0.09 


-9.05 ±0.05 


-8.81 ±0.05 


-8.93 ±0.03 


-8.77 ±0.08 


-8.85 ±0.10 


-8.55 ±0.09 



"Empirically measured per pixel at 3250A. 
^Units of ltr 5 
c Weighted mean 

References. — Key to References - 1: Laurent, Vidal-Madjar, & York (1979); 2: Jenkins et al. (1999); 3: Hoopes et al. (2003); 4: 
Sonneborn et al. (2002); 5: Friedman et al. (2002) 



In this paper, we report on our first set of Ti II obser- 
vations of seven sightlines. We measure the Ti + column 
densities to assess the depletion levels along the sightlines 
and examine correlations with the observed D/H values. 
Finally, this Letter establishes a public database for Ti II 
measurements obtained by our group 5 . The data pre- 
sented here, and all future observations, will be archived 
at this site, including the raw data and calibration frames. 
The data are freely available to other researchers. 

2. OBSERVATIONS AND REDUCTION 

The seven sightlines presented here were observed on 
the nights of September 8 and 9 and October 6, 2004 UT 
with the recently upgraded HIRES spectrometer (Vogt et 
al. 1994) on the 10m Keck I telescope. The spectrome- 
ter now contains a three CCD mosaic with > 90% quan- 
tum efficiency at A < 3500A. Because the targets are 
bright, we obtained the data primarily during twilight. 
The first half of the first night was marred by poor ob- 
serving conditions while the remainder of the time was 
clear with typical seeing of FWHM k, 0.7". The targets 
in September 2004 were observed through the B5 decker 
(0.86" width; FWHM w 6km s" 1 resolution; 3.5" slit 
length) while BD±28°4211 was observed through the Bl 
decker (0.57" width; FWHM « 4.5km s" 1 resolution; 3.5" 
slit length). Table 1 summarizes the exposure times and 
data quality of the sample. 

At the time of publication, a data reduction pipeline for 
the upgraded HIRES instrument was not available. There- 
fore, we reduced the two echelle orders containing the Ti II 
A3230, 3242, 3384 transitions with a set of in-housc IDL 
routines. These routines subtracted the bias, interactively 
set a boxcar aperture, interactively determined a region for 
scattered light subtraction (sky background was minimal) , 

5 http://www.ucolick.org/-~xavier/TiII/index.html 



and extracted a ID spectrum. A variance array was calcu- 
lated accounting for the read noise and assuming Gaussian 
statistics. Wavelength calibration was carried out by fit- 
ting a 3rd order polynomial to the arc lines identified in 
the ID spectra extracted from a ThAr image using the 
same trace and aperture as the science extraction. The 
typical RMS deviation from the fit was less than 0.1 pix- 
els (i.e. < 0.003A). Although we achieved a formal S/N 
ratio > 500 pix -1 for the brightest targets, the scatter 
around the stellar continuum significantly exceeds Gaus- 
sian statistics in the brightest sources. This is most likely 
due to noise in the combined flat field or perhaps small er- 
rors in the traces of the dispersed spectra. Throughout the 
analysis, we have augmented the variance arrays to match 
the empirical scatter measured in the stellar continuum. 

All of the individual exposures were rebinned to a com- 
mon wavelength scale with A = 3000A and a pixel size 
of 1.4km s _1 after correcting to vacuum wavelengths and 
the heliocentric velocity reference frame. Multiple expo- 
sures of the same object were compared to identify cos- 
mic rays and then coadded after scaling to a common flux 
and weighting by the square of the median signal-to-noise. 
Finally, we normalized the data by fitting a high order 
polynomial to the unfluxed stellar continuum. Continuum 
normalization is generally the largest source of uncertainty 
in the analysis, especially for the weakest transitions. 

3. ANALYSIS 

Figure 1 presents a velocity plot for the HD 191877 
sightline. At the recorded S/N and resolution, it is straight- 
forward to measure the Ti II ionic column densities. In 
the following, we assume that the majority of Ti+ is in 
the 4 F 3 / 2 ground state. From our observations, we mea- 
sure A^(TiII*A3239) < 10 10 - 6 cm~ 2 as a conservative up- 
per limit to the 4 F 5 / 2 state (E « 135 K). Provided the 
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Fig. 1.— Till profiles for the ISM sightline toward HD 191877. 



Ti II profiles have small optical depth, we can measure the 
Ti II column density by (i) summing the EW and assume 
the weak limit for the curve-of-growth (the values range 
from several to 30 mA) ; (ii) integrating the line profile with 
the apparent optical depth method (AODM; Savage and 
Sembach 1991; Jenkins 1996); and (iii) by fitting Voigt 
profiles to the data. All three techniques yield similar re- 
sults; we have mainly used the AODM in this paper. The 
results for the three transitions 6 for all of the sightlines 
are listed in Table 1. For all measurements we report la 
uncertainties and the upper limits correspond to 3cr lim- 
its. The uncertainties include (in quadrature) statistical 
uncertainty from Poisson noise and error due to contin- 
uum placement. We estimate a k continuum error of 
0.1% for all of the objects except FeigellO (0.2%) which 
corresponds to an approximately 0.5 mA error over the in- 
tegrated profiles. Finally, we adopt a minimum error of 
0.02 dex owing to systematic uncertainties related to data 
reduction (e.g. flat fielding). 

Figure 2 presents log(Ti/H) against the (D/H) values 
reported in the literature. The visual impression is sug- 
gestive of a correlation between the two quantities. The 
non-parametric Spearman and Kendall correlation tests 
reject the null hypothesis at the 96% and 95% c.l. respec- 
tively and indicate a positive correlation between the two 
quantities. This correlation provides preliminary evidence 
that the variation in Galactic D/H is related to the phys- 
ical conditions that imply high depletion levels. Indeed, 
the results are consistent with the interpretation that low 
Galactic D/H values are due to an increased depletion of 
D onto grains. 

Examining the figure, it is obvious that additional mea- 

6 Assuming / 323 o = 0.0687, / 32 42 = 0.232, / 3384 = 0.358 (Bizzarri et 
al. 1993; Pickering, Thorne, & Perez 2001, 2002; Morton 2003). 
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Fig. 2. — log(Ti/H) vs. D/H for the seven sightlines comprising 
the current sample. Error bars reflect la uncertainty. 



surements of sightlines with large (D/H) value will be 
especially valuable for testing the correlation. Further- 
more, one observes a large spread in Ti/H at low D/H. 
If confirmed by future observations, the spread in Ti/H 
would argue against a simple linear correlation between 
Ti/H and D/H. We note that two additional stars (£Pup, 
7 2 Vcl) with (D/H) > 1 x 10 -5 have Ti II column densi- 
ties reported in the literature (Welsh et al. 1997). The 
inferred Ti/H value for (Pup (log(Ti/H) = -8.5 ± 0.08) 
follows the trend indicated in Figure 2. The upper limit to 
Ti/H (< —8.9) for j 2 Vel, however, suggests the gas along 
this sightline is significantly depleted even though its D /H 
value is among the highest known. A similar conclusion 
may be drawn from the GHRS observations of Fe II by 
Fitzpatrick & Spitzer (1994) although line-saturation is 
a potential concern. Although the confirmation of a high 
depletion level toward 7 2 Vel would raise concern, Draine 
(2004a) notes that the processes of D depletion are differ- 
ent from those for other ions, e.g., D may be depleted in 
polycyclic aromatic hydrocarbons while heavier elements 
are depleted by other types of dust grains. Therefore 
one need not expect a one-to-one correspondence between 
Ti/H and D/H. In the coming year, we plan to acquire 
Ti II observations along many additional sightlines to bet- 
ter constrain the slope and scatter of Ti/H vs. D/H. 

We have also examined the correlation between Ti/H 
and the iV(HI) values of the sightlines. Because one ob- 
serves a correlation between volume density and deple- 
tion level (e.g. Jenkins 1987), one may expect a similar 
trend for A r (HI). Furthermore, if photoionization is impor- 
tant along these sightlines its effects should correlate with 
AT (HI). The Spearman and Kendall tests, however, reveal 
the null hypothesis is ruled out at only the 29% and 45% 
c.l. respectively. This follows the results from previous 
surveys of the Galactic disk Ti II (Stokes 1978; Welsh et 
al. 1997) although Wakker & Mathis (2000) see a strong 
correlation in a set of clouds with much larger dynamic 
range in iV(HI) than the sightlines considered here. 

We wish to emphasize an important aspect of differen- 
tial depletion in the analysis of Galactic D/H. In Figure 3, 
we present the Ti II profiles for the HD 195965 sightline 
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Fig. 3.— Velocity profiles for the ISM sightline toward HD 195965. 
Separately, the refractory ions (Ni+, Ti+) and non-refractory ions 
(P + , O ) show very similar sets of profiles. One notes, however, 
that the positive velocity component has a much larger depletion 
level for Ni+ and Ti+ (m 10 x) than the negative component. 



against a series of UV transitions obtained with STIS on 
the Hubble Space Telescope. The UV data have higher 
resolution (FWHM w 1.5km s _1 ) but lower S/N per res- 
olution element than the optical observations. Examin- 
ing the detailed component structure of the profiles, the 
refractory ions (e.g. Ti + , Ni + ) have similar characteris- 
tics. Comparing against the non-refractory profiles (P + , 
O ), however, we identify the two main components but 
note that the relative abundances are significantly differ- 
ent. The Ti abundance relative to O is 10 x greater in 
the positive component than the negative component. It 
is also noteworthy that the Mg II profiles more closely 
track the non-refractory species even though Mg is de- 
pleted along the sightline. 

This raises a number of concerns regarding the analy- 
sis and interpretation of D/H. First, variations in the de- 
pletion level on a given sightline would generally lead to 
a weaker correlation between any existing Ti/H vs. D/H 
correlation if one only considers the integrated values. By 
a similar token, one would tend to underestimate the mag- 
nitude of intrinsic scatter in D/H regardless of the physical 
mechanism responsible. Another complication is that it is 
unclear whether one should constrain the analysis of D in 
lower resolution data (e.g. FUSE observations) with the 
velocity profiles traced by the O I, Mg II, or Ti II profiles. 
Profiles of Fe II are likely to be similar to Ni II and Ti II 
and could also be confusing in analysis of D I. A useful ex- 



ercise would be to quantitatively compare high resolution 
D profiles with both refractory and non-refractory transi- 
tions along sightlines where the depletion levels vary. We 
consider these issues and examine the systematic effects 
on the inferred D/H ratios in a future paper. Finally, we 
note that the D profile of at least one extragalactic sight- 
line shows significantly different component structure than 
the corresponding metal-line profiles (Tytler et al. 1996). 
Although the physical origin is more likely related to dif- 
ferences in metallicity or the ionization state of the metals, 
one cannot rule out intrinsic scatter in D/H even in these 
low metallicity sightlines. This issue could be relevant in 
studies of D/H in damped Lya systems. 
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NSF grant AST-0307408. TMT appreciates support from 
NASA through grant NNG04GG73G. This research was 
supported in part by the National Science Foundation un- 
der Grant No. PHY99-07949. 
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